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Abstract 

We consider the viability of dark energy (DE) models in the framework of 
the scalar-tensor theory of gravity, including the possibility to have a phantom 
DE at small redshifts z as admitted by supernova luminosity-distance data. 
For small z, the generic solution for these models is constructed in the form of 
a power series in z without any approximation. Necessary constraints for DE 
to be phantom today and to cross the phantom divide line "p = —p at small z 
are presented. Considering the Solar System constraints, we find for the post- 
Newtonian parameters that 7p7v < 1 and ^pn,o ~ 1 for the model to be viable, 
and l3pNfi > 1 (but very close to 1) if the model has a significantly phantom DE 
today. However, prospects to establish the phantom behaviour of DE are much 
better with cosmological data than with Solar System experiments. Earlier 
obtained results for a A-dominated universe with the vanishing scalar field 
potential are extended to a more general DE equation of state confirming that 
the cosmological evolution of these models rule them out. Models of currently 
fantom DE which are viable for small z can be easily constructed with a constant 
potential; however, they generically become singular at some higher z. With a 
growing potential, viable models exist up to an arbitrary high redshift. 
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1 Introduction 



A major turning point in cosmology has been reached with the observational discovery 
that our Universe is accelerating now (and has been accelerating for several billion 
years in the past) [HE]- If interpreted in terms of the Einstein equations for the 
evolution of a Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmological models 
with the (practically) zero spatial curvature (the latter follows from other arguments), 
this means that approximately two thirds of the total present energy density of matter 
in our Universe is due to some gravitationally unclustered component called Dark 
Energy (DE). 

Observations of high redshift supernova, fluctuations of the cosmic microwave 
background (CMB) temperature and other effects tell us that the effective energy 
density poE of DE is very close to minus its effective pressure pde (see Eqs. f lTGlTTl) 
below for the exact, though conventional, definition of these quantities) and that they 
are both very weakly changing (if at all) with time and with the expansion of the 
Universe. The physical nature of DE is unknown at present, and three main logical 
possibihties exist (see [21 HJ [5l El [7] for reviews). 

1. DE is a cosmological constant, as originally suggested by Einstein, with poE = 



2. Physical DE: DE is the energy density of some new, very weakly interacting phys- 
ical field (e.g., a quintessence - a scalar field 4> with some potential 1^(0) minimally 
coupled to gravity). 

3. Geometric DE: the Einstein general relativity (GR) equations are not the correct 
ones for gravity, but we write them in the Einsteinian form by convention, putting all 
arising additional terms into the r.h.s. of the equations and calling them the effective 
energy-momentum tensor of DE. 

Of course, this classification is not absolute. In some cases, the difference be- 
tween physical and geometric DE can, in turn, become conventional. E.g., for a 
non-minimally coupled scalar field (which constitutes a specific case of scalar-tensor 
gravity considered in this paper), equations of the model have the same form irre- 
spective of the origin (physical or geometric) of this field. Another way of classifying 
DE models which is more invariant and, therefore, more important from the point of 
view of the observational determination of the type of DE is to divide these models 
into dynamical, if the DE description requires a new field degree of freedom (= a new 
particle from the quantum point of view), and non-dynamical in the opposite case. 
Then all physical DE models and many of the geometric ones belong to the first class, 
while the cosmological constant itself and some geometric models fall into the second 
category, i.e., the F{R) model with the Palatini variation of its action where R is the 
scalar curvature. 

At present, all existing observational data are in agreement with the simplest, 
first possibility (inside ~ 2a error bars in the worst case). The case of a cosmological 
constant is internally self-consistent and non-contradictive. The extreme smallness of 
the cosmological constant expressed in the either Planck, or even atomic units means 

= c = 1 is used throughout the paper. 
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only that its origin is not related to strong, electromagnetic and weak interactions (in 
particular, to the problem of the energy density of their vacuum fluctuations) at all. 
However, in this case we remain with one dimensionless constant only and can not say 
anything more (at least, at present). That is why the two other possibilities admitting 
a (slightly) variable dark energy have been also actively studied and compared with 
observational data recently. Moreover, properties of the present DE are remarkably 
qualitatively similar to those of an "early DE" that supported an inflationary stage 
in the early Universe. But in the latter case, we are sure that this early DE was 
unstable. So, it is natural to conjecture by analogy that the present DE is not stable, 
too. 

On purely phenomenological grounds, one can consider DE models with a constant 
equation of state parameter wde = Pde/pde different from —1. However, the latest 
observational data have already severely restricted this simplest alternative possibihty 
to + ~ 0.1 (Icr error bars) [Bl El [ID]- Therefore, a viable alternative to the 
cosmological constant has to be looked for among more complicated models with 
wde ^ const. Such models include quintessence ones with different potentials V{(f)) 
(see [11 mini [6] for numerous references), models with several minimally coupled scalar 
fields (see, e.g., [TTj), those with direct non-gravitational coupling between DE and 
dark matter (DM) ([12] and following papers), unified models of DE and DM (the 
Chaplygin gas model [T3] and others), etc. Then, however, it becomes very important 
to investigate if there exist models of DE for which a variable wde may cross the 
"phantom divide" line wde = —1 (we call DE "fantom" if it has wde{z) < —1 for 
a given z and "normal" in the opposite case). Note that the weak energy condition 
(WEC) is violated for phantom DE. This is not possible at all for quintessence models 
with the standard kinetic term and hardly possible, i.e., it requires non-generic initial 
conditions, for scalar field models having a non-standard kinetic term (" /c-essence" ) 
[m (see also [15]). 

Indeed, analysis of the recent SNe data with redshifts up to 2; = 1.7 (the "Gold 
dataset" [IB]) using fits containing at least 2 free parameters, e.g., the linear fit 
( 175]) for Wde in terms of the scale factor a [ITJ dS] or the quadratic polynomial fit 
f l77|) for pde in terms of the redshift z [19], or with model independent methods, 
results in best fits to these data having a variable wde that steadily increases for 
redshifts < z < 1 and crosses the "phantom divide" somewhere between and 
0.5 [SniEIlElESlEaESlESlETlESlEn] (see, however, [SUl El EE E2] for a more 
conservative view - in the sense of returning back to a cosmological constant). This 
statement does not mean that an exact cosmological constant is excluded - it is still 
inside 1 or 2a error bars, as was stated above. Moreover, one should be cautious 
with this result: it may be a consequence of trying to obtain a too fine-grained graph 
of wde{z) from the luminosity distance dii^z) as has been already emphasized in 
[53] . In this respect, results for wde{z) averaged over a range of redshifts > 0.3 
(dubbed the "w-prohe" in ^3]) are more statistically reliable, e.g., the result that 
Wde ~ — 1-05]'^o;q7 obtained for the same Gold SNe dataset in pm for fi^.o = 0.3 
and the averaging redshift interval (0,0.414). However, a remarkable possibility of 
crossing the phantom divide at recent redshifts, when DE has become the dominant 
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component of the Universe by its effective energy density, still remains viable. The 
latest Supernovae data with z ~ 0.5 [35] also admit wde{z) < —1 for z < 0.5, but not 
for larger redshifts, as a possible interpretation. Finally, the recent data on acoustic 
baryon oscillations in the present matter power spectrum (the Sakharov oscillations) 
[36l [37] . while strongly restricting one direction in the plane of parameters for the 
two-parametric fits fl76f77p . leave the orthogonal direction practically unconstrained, 
therefore, permitting the recent phantom divide crossing (see Sec. 6). 

So, if this striking behaviour of DE will be confirmed by future data, what is the 
best way to describe it? One possibility, a ghost phantom DE, i.e. a scalar field with 
the negative sign of its kinetic term, was first proposed in [38] and triggered a very 
large wave of publications (see, e.g., the recent papers [321 SO] for a list of references 
on this topic). However, it has been long known that theories of this type are plagued 
by quantum instabilities, the most dangerous of those being the process of creation of 
two particle+antiparticle pairs: one of the ghost field and another of any usual (non- 
ghost) field (see [31] for a recent investigation). Moreover, a ghost model of DE is 
unsatisfactory even at the classical level: it does not explain the observed large-scale 
isotropy and homogeneity of the Universe! Just the opposite, e.g., for the given Hubble 
constant Ho averaged over angular directions in the sky, we would expect a universe 
to be very strongly anisotropic with the anisotropy energy density (i.e., the positive 
energy density of long-wavelength gravitational waves) being compensated by a large 
negative energy density of the ghost DE, or of a ghost component of DE in more 
complicated multi-component models of this type (e.g., in the two-field realization 
of the "quintom" model introduced in [12]). It is just a classical analogue of the 
quantum instability mentioned above, with the "usual" field being the gravitational 
one and with the dynamical quantum instability transformed into the problem of 
classical initial conditions. For this reason, we are sceptical regarding this approach 
whole. 

Fortunately, there is no necessity to introduce a ghost field to explain possible 
phantom behaviour of DE including its phantom divide crossing. As was first empha- 
sized in [33] , scalar-tensor theories of gravity allow for this phenomenon. Scalar-tensor 
models of DE belong to the third class (geometric DE) according to the classifica- 
tion given above. This class of models is very rich and interesting. It contains the 
Einstein gravity plus a non-minimally coupled scalar field, as well as the higher- 
derivative f{R) gravity, where R is the scalar curvature, as particular cases, besides 
allowing for DE phantom behaviour and transition to a normal DE. In this paper, 
we will concentrate on scalar-tensor DE models with the latter properties because 
of tentative observational evidence described abovej^ Note that a possible phantom 
behaviour of DE in scalar-tensor gravity has a conventional character. The reason is 
that the effective gravitational constants Gn and Gefr (see Sec. 2 below) are generi- 
cally time-dependent while the definitions fll6|17p of the DE effective energy density 
and pressure assume some constant G when writing the left-hand side of equations 
and splitting their right-hand side into energy-momentum tensors of non-relativistic 

■^The third possibility to get a phantom DE including the phantom divide crossing which is based 
on braneworld models [44] (see also the recent paper [45]) will not be discussed here. 
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matter (mainly non-baryonic cold dark matter) and dark energy. As a result, a part 
of the Einstein tensor G^^ and the energy-momentum tensor of dark matter multi- 
plied by a change in Gat is conventionally attributed to DE. In other words, phantom 
behaviour of scalar-tensor DE is always 'curvature- induced', in contrast to the ghost 
DE models or other models like those considered in [l6l HTj where it may occur in 
the flat space-time already. 

The possibility to get both a phantom DE and the phantom divide crossing in 
scalar-tensor gravity is related to the fact that this theory has two arbitrary and inde- 
pendent functions and f/($) (see the Lagrangian ([1]) below). Throughout this 
work we assume spatial flatness though this prior while well motivated theoretically 
can be challenged [48j. As has been shown in ^43j, two different types of observations, 
e.g., determination of the luminosity distance and the inhomogeneity growth factor 
in the non-relativistic component as functions of redshift, are necessary and sufficient 
for the total reconstruction of the microscopic Lagrangian ([1]) of scalar-tensor gravity. 
However, as shown in po], a partial reconstruction using the luminosity distance data 
only could yield interesting information, too. In that case some assumption about 
one of the functions F{^) and t/($) has to be made, so that only one unknown func- 
tion in the microscopic Lagrangian ([1]) remains to be found. Clearly, many different 
partial reconstruction strategies are possible and we explore some of them here with 
the aim to investigate which DE models in scalar-tensor gravity are viable. 

It has been found in flS] that A-dominated universes with a vanishing potential U 
are ruled out as they lead to singular universes already at very low redshifts z ~ 0.7 
(see also the recent paper [50]). Of course, completely regular but non-accelerating 
solutions do exist in this case. Among models with a non-zero potential U, one 
concrete example of a model with the phantom divide crossing was constructed in 
[50\ where the functions -F($) and U ($) were given in the parametric form as functions 
of z up to 2 ~ 3. More examples for a non-minimally coupled scalar field {F{^) = 
Fq — ^$^) with a non-zero potential U were investigated in [SI], also for z < 2. In 
the present paper, we make a next step and extend these results in two directions: 
first, by constructing a generic solution for scalar-tensor DE models for z < 1 in the 
form of a power series in z; second, by investigating and numerically integrating some 
of these solutions up to large z ^ 1. The latter task appears to be necessary since, 
though DE is subdominant for z ^ 1 as compared to non-relativistic non-baryonic 
dark matter and baryons, its model itself may become intrinsically contradictory for 
large z, namely, F or $^ may become negative for an unfortunate choice of f/($). It 
is also crucial to check that any DE model has the correct power-law behaviour for 
large z [52] . 

Finally, for the scalar-tensor theory of gravity, it is well known that the present 
value of ^(2; = 0) is severely restricted from Solar System tests of post-Newtonian 
gravity (i.e. by the measured values of the post-Newtonian parameters). For this 
reason, there have been stated that the present phantom behaviour in scalar-tensor 
models of DE requires large amount of fine tuning and, thus, is unnatural [53l 
Therefore, it is important to investigate this problem in more detail to quantify what 
amount of fine tuning (and of what kind) is necessary for a significantly phantom 
behaviour of DE, and to determine the relation between this behaviour and results 
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of Solar System tests of gravity. 

The paper is organized as follows. In Sec. 2 we define all quantities related to our 
scalar-tensor DE model and present the background evolution equations. 

In Sec. 3 we derive the general integral solution for the quantity F{z). In Sec. 4 
we consider solutions in which DE scales as some power of the FLRW scale factor a{t) 
and show their existence for DE of the phantom type (the latter requires a non-zero 

In Sec. 5 we consider the small z behaviour of our model and find conditions for the 
violation of the WEC by DE today and for the phantom boundary crossing at small 
z. In Sec. 6 the general reconstruction of a microscopic model is considered and the 
observational constraint from the acoustic oscillations in the matter power spectrum 
is discussed. In Sec. 7 we consider the reconstruction for a constant potential more 
specifically and show that the model becomes singular at some redshift that cannot 
be arbitrarily high once we have chosen a specific equation of state wde- In Sec. 8 
non-constant potentials are considered and a model which is asymptotically stable 
for large z is presented. Sec. 9 contains conclusions and discussion. 



2 Background evolution 

In this section we review the background evolution equations in a spatially flat FLWR 
universe. We consider a model where gravity is described by a scalar-tensor theory 
and we start with the following microscopic Lagrangian density in the Jordan frame 

L = i R - Z($) g^''d^<!>d,<^^ - f/($) + Lm{g,u) • (1) 

Since Lm is not coupled to the Jordan frame is the physical one. In particular, 
fermion masses are constant and atomic clocks measure the proper time t in it. The 
quantity Z{^) can be set to either 1 or —1 by a redefinition of the field $, apart 
from the exceptional case Z{^) = when the scalar-tensor theory ([T]) reduces to 
the higher-derivative gravity theory R + f{R)- In the following, we will write all 
equations and quantities for the case Z = 1. For our purposes, Lm describes non- 
relativistic dust-like matter (baryons and cold dark matter) as we are interested in low 
redshift {z -C Zeq) behaviour only. Here, Zeq denotes the equality redshift when the 
energy densities of non-relativistic matter and radiation are equal. In such a model, 
the effective Newton gravitational constant for homogeneous cosmological models is 
given by 

= i8nF)-\ (2) 

As could be expected, does not have the same physical meaning as in General 
Relativity, the effective gravitational constant Gefr for the attraction between two test 
masses is given by 



F + 2{dF/d^) 



2 



on all scales for which the field $ is effectively massless [13] and F > 0. The condition 
Geff > is one of the stability conditions of the scalar-tensor theory ([1]), it means 
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that the graviton is not a ghost. In fact, even at the purely classical level, it has been 
shown in [55] that a generic solution of may not smoothly cross the boundary 
Ges = 0. Instead, a curvature singularity forms at this boundary which generic 
structure has been also constructed in [55]. This condition combined with another 
stability condition (see Eq. (fTOj) below) results in F > 0, so > 0, too. 

We will write all equations in the Jordan frame using ([T]). Specializing to a spa- 
tially flat FLRW universe 



-df + a^(t) dx^ 



the background equations are as follows: 



3FH^ 
-2FH 



Pm 
Pm 



$2 

+ — + U- 3HF 
+ + F - HF . 



(4) 

(5) 
(6) 



The evolution equation of the scalar field $ can be obtained from the two equations 
( l5f6l) . Eliminating the quantity by combining these equations, we obtain a master 
equation for the quantity F which takes the following form when all quantities are 
expressed as functions of redshift z: 



F" + 



(In/i)' - 



l + z 



F' + 



6 



'l + z) 



2{\nhy 

l + z 



Qu 



[l + zY 



Fo nu,o + 3 (l + z) h-^Fo a 



in,0 5 



(7) 



where a prime denotes the derivative with respect to z and we have introduced the 
quantities h = ilu^ = g^,^"^^ and u = The index denotes the present 

moment here and below. The (dimensionless) relative energy density fim.o is defined 
through flmfi = 3^^- Once the master equation ([7]) is solved for F, we get the 
algebraic equation for $'(2): 



$'2 



F' 



+ 



Fn u 



Fo (1 + z) 



l + z ' (1 + ^)2 (1 + Z)2 ^'^ 

The second stability condition of the scalar-tensor gravity ([T]) is 



m,0 



/2 



^BD 



p/2 



> 



3 
2 



(9) 



where ^'^ is found from ([8]). Inequality Qj just expresses the positivity of the energy 
of the (helicity zero) scalar partner of the graviton, i.e. the positivity of the kinetic 
energy of the scalar field in the Einstein frame (see e.g. [49] for more details) 



A'2 



F' 



$'2 

— > 
2F 



(10) 
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with taken from (jH]). With the Z = \ parametrization, we cannot reconstruct the 
function when — | < ujbd < since becomes negative in this case. Indeed, 

these allowed negative values ofusD correspond to the parametrization choice Z = —1 
in Eq.([T]). The Z = 1 parametrization allows us to consider consistently only cases 
for which > 0, or equivalently ujbd > 0. However, the condition (Q with 
given by ([H]) is the true condition that the theory is well behaved and it remains 
valid even for < 0. This is best understood in the Brans-Dicke parametrization 
F = $, Z = , where one can reconstruct the two functions F > and ujbd > — | 
from ([7]), ([8]), 0. For > and F > 0, the inequality is satisfied automatically. 
The inclusion of the range —^<ubd<^^s not of purely academic interest and can 
be important when one considers the reconstruction of DE models far in the past. 
We illustrate this with a specific example in Fig. 1. 

Solar System experiments constrain the post-Newtonian parameters ■jpn and Ppn 
today (for these quantities, we drop here the subscript 0) 



F + 2{dF/d^Y ' ^ ' 

_ 1 F (dF/d^) rf7 



G 

as well as the quantity The best present bounds are: 

= (2.1 ±2.3) ■ 10"^ 
l3pN-l = (Oil) -10-" 

^ = (-0.2 ±0.5)- 10-13 (13) 
<-^efr,o 

where the first bound was obtained from the Cassini mission [56] and the other two 
from high precision ephemerides of planets [57] (the second bound has been recently 
confirmed by the Lunar Laser ranging [58] -their value is/3pAr — 1 = (1.2±1.1) ■ 10~^). 

As a consequence of the smallness of jpn — 1, the Brans-Dicke parameter ubd 
satisfies today the inequality 

iUBDfi > 4 X 10^ (14) 
The resulting bound on F'{0) is very stringent 



F'(0) ^ $'(0) 



(15) 



However, as was discussed in [13], the quantity ojbd need not be so large as (|T4l) in 

the past, though one can deduce a looser inequality applying up to redshift ^ ^ 1 

with fairly reasonable assumptions. From ([8], [T5l) we can derive the allowed range of 

_ i_ 

initial values F'(0) and we find |F'(0)| ~ tu^l, q and ^oEfi — ^ufl > 0, the result that 
we will recover in Sec. 5 when performing the small z expansion of all quantities. The 
peculiar case F'(0) = together with Qde,o ~ ^u,o > corresponds to pure General 
Relativity today {ujbd = oo)- 
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As noted above, supernova observations permit DE to be of the phantom type, 
with the equation of state parameter wde < —1 at small redshifts. This is a strong 
motivation for considering DE models in the framework of the scalar-tensor theory 
of gravity. More generally, at present there is much interest in models with modified 
gravity of which the scalar-tensor theory is a well known representative. In scalar- 
tensor DE models, a meaningful definition of energy density and pressure of the DE 
sector requires some care (see also [,53j for a detailed explanation). Let us define the 
energy density pde and the pressure pde the following way: 

3FoH^ = Pm + PDE (16) 

-2Fq H = pra + Pde + Pde ■ (17) 

This just corresponds to the (conventional) representation of the true equation for 
scalar-tensor gravity interacting with matter in the Einsteinian form with the constant 
Co = GNito)'- 

Rfiu — 2^3fiu = SttCq iT^u,m + T^u,de) ■ (18) 
With these definitions, the usual conservation equation applies: 

Pde = -3H{pde + Pde) ■ (19) 
If we define the equation of state parameter wde through 

Wde = , (20) 

Pde 



the time evolution of the DE sector is given by 

pDEfl I Jo I + Z' 

Using fll6|17l) . one gets 
where 



(21) 



2/1 I ^\dlnH 



:i+^) 



WDE = , (22) 

1 - ^^,n,o(l + ZY 
^" ^ • ^^^^ 

Then Eq.( ll6p can be rewritten as 

h\z) = fi^,o (1 + Zf + ^DEfi <z) (24) 

where ^dea = 1 ~ ^m.o by definition. The condition for DE to be of the phantom 
type, Wde < —1, can be obtained from (!22|) . It reads [3l H3] 

^<3n^{l + zf . (25) 
dz 
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This inequality is modified in the presence of spatial curvature |1S]. As was first 
emphasized in [13] , the scalar-tensor gravity allow for phantom DE. Indeed, for these 
models 

PDE+PDE = ^^ + F-HF + 2{F-Fo) H , (26) 

hence the weak energy condition for DE can be violated (see also [53]). Moreover, 
as it will be shown below, the weak energy condition may be violated even for the 
sum of DE and non-relativistic matter, i.e. for the whole right-hand side of Eqs. 
( ITSjl . leading to dh/dz < 0. However, such a strong violation corresponding to wq = 
wde{0) < (1 — Qmfl)^^ ~ —1.4 is not supported by the existing data (though not 
completely excluded either). 

The relation between the Hubble parameter and the luminosity distance in scalar- 
tensor gravity is the same as in GRq 

h''(') = {j^) ' d^iz) = Ho Dl{z) . (27) 

However, as discussed in [59], [60], when using supernova data to obtain Dl{z), one 
has to take into account the dependence of the Chandrasekhar mass on Gcfr, so that 
the Snia peak luminosity appears to be oc (G'efr(-2)/G'cfr,o)~^^^- As shown in |;43j, the 
full reconstruction of the functions -F($), U{^) requires two independent types of 
observations: dL{z) or another function of z that probes the background evolution, 
and the growth factor of matter perturbations Spm{z)/pm{z) at some comoving scale 
much less than the Hubble scale, too. On the other hand, as emphasized in [19], one 
can already obtain powerful constraints from a partial reconstruction using Dl{z) 
only. Such a partial reconstruction is possible when some additional condition is 
imposed on either F or U. This is the way we adopt in Sec. 7: we reconstruct 
the function F for a given Hubble parameter H{z) and a constant potential U and 
investigate whether the resulting model is viable. Before embarking on such partial 
reconstructions, we derive first an integral form of the general solution of Eq. ([7j). 



3 The master equation for F 

Let us now consider the master equation ([7]) and present its general solution in an 
integral form. This will allow us to understand general properties of solutions, in 
particular, their dependence on initial conditions. The first step is to note that one 
solution of the homogeneous equation (iTl) (without source term) is given by 

Foc{l + zf. (28) 

This suggests us, following [49], to introduce the function / defined as follows 

^P^(l + z)V(l + ^), (29) 



■^Note that the definitions of dL{z) and Dl{z) are interchanged in [T7] and [15] . 
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in terms of which Eq.d?]) becomes 



r + {\nhyf 



6 u 



X- 



^U,0 + 



X h"^ 



m,0 ; 



(30) 



where we have introduced the variable x = 1 + z. Due to the absence of any term 
proportional to /, in accordance with fl28|29p . Eq. (1301) is easily integrated formally, 
and its general solution has the integral form 



fix) = 1 + 



-2 + 



F'(x = 11 



~ ^ ii 7^0^ ii x"4 h{x") 
r dx' r' dx" 



A nontrivial dependence on initial conditions is through the second term only. 
We see that for a given initial value F'{z = 0), this term is proportional to the 
dimensionless luminosity distance d^. Due to Eq. (1151) . we have f'{x = 1) = —2 ± 
\j> (2-0)1 ^ ^ indicating that we must start today, on observational grounds, very 

close to GR. Hence, the second term of the general solution is bound to be negative 
and it is this term that will possibly push the quantity f{z) downwards for increasing 
z. Finally, the corresponding quantity F is trivially obtained using (l29l) . If we choose 
U to be constant, we can implement the reconstruction of the microscopic Lagrangian. 
By inspection of ([7]) for constant U, we find the asymptotic solution 



F{x) =CiX^ + C2 x-2 + Fo 



(32) 



in complete agreement with (|3T|) setting u = 
reconstruction scheme |l3l HQ] and find $(z 
This would finally give us F{^). 

Let us come back to the solution 
Z = —1 with 



1. We can now proceed with the general 
, and hence z{^), by integration of ([8]). 

It corresponds to the parametrization 



3 

"2 



(33) 



where $ is defined up to some constant. From (1331) it is seen that this solution 



is unphysical in view of (19]), see [19]. It is interesting that (13311 corresponds to a 
conformally coupled scalar field in the Jordan frame satisfying 



$ + 3H^ + -$ 
6 



. 



(34) 



4 Scaling solutions 



Let us now consider the so-called scaling solutions for which DE scales as some power 
of a, 

Pde oc a~^^ , (35) 
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which attracted a lot of interest previously!^ Clearly, for these solutions DE has an 
effective barotropic (constant) equation of state, viz. 

wde = -1 + 7 , (36) 
and we can readily write formally the general scaling solution by substituting 

= nm,o (1 + + nnEfl (1 + zf^ , (37) 

into the integral expression (l3Ti) . We emphasize that this gives the general scaling 
solution irrespective of any limiting case and for an arbitrary potential, the only 
assumption is that of the spatial flatness. Using (!3T|) and (137|) . general analytic 
expressions for F{z) can be obtained only for the cases 7 = 1, 4/3 (see Appendix). 

In analogy with a minimally coupled scalar field (quintessence), a scaling solution 
satisfies 

Pde + Pde ^ , .„„^ 

= 1 + Wde = 7 • (38) 

Pde 

However, in contrast to the minimally coupled scalar field case with a positive poten- 
tial, (|38|) does not imply < 7 < 2, in particular 7 can be negative which corresponds 
to phantom DE. 

To get insight into the ability of scalar-tensor DE models to produce various 
equation of state parameters wde, it is instructive to study first scaling solutions 
in the absence of dust-like matter (p^ = 0). In this way, a lower limit on wde for 
realistic solutions with Pm 7^ can be obtained. These scaling solutions can also be 
considered as describing the asymptotic future of our universe when and 

^Ide — > 1- Then a{t) oc \t\'^ with q = For phantom (or, super-inflationary) 
solutions, g < and then the moment t = corresponds to the 'Big Rip' singularity 

[all Eg. 

It is straightforward to check that such scaling solutions can exist only when 
F = 0;$^ with a = const. This form of F corresponds to a constant Brans-Dicke 
parameter ujbd = 4^ for 2'($) = 1. Further on, we assume that a > for stability 
of the theory. The two scaling solutions in the vanishing potential case (the pure 
Brans-Dicke theory) are 



l±3Jl + ^ l + J-^Jl + ± 

Hence it is not possible to get a scaling solution with woe < —1 with vanishing 
(or negligible) potential U{^). We note that the above solution reduces to (I28II33P 
obtained for Z = — 1 when we formally put a = — | in (15^ . Indeed we get q = I, r = 
— 1, hence F oc which is exactly the solution (1281) . 

However, scaling solutions supported by both the kinetic and the potential energy 
of $ can exist in the presence of a polynomial potential f/($) = f/o|$|" (we assume 



^Our definition for pde, Eq. p6p . differs from that used in [61j . and this explains why our results 
differ from theirs. 
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that the potential is positive, so Uq > 0). We have for these solutions 



_2 _2(n + 2 + i) 

~ 2^ ' ^ ~ {n-2){n-4) 

while the inequality + 2qr > ^ has to be satisfied (this is always the case if a ^ 1 
and n = 0{1)). They were first found in [Bl] (see also the analysis of their stability 
in ES]) but they have not been discussed yet in connection with phantom DE. 
Note that there is no phantom behaviour for usually considered quadratic and quartic 
potentials. 

So, we see that it is possible to get scaling solutions with g<0, wde < 
if 2 < < 4. For these solutions, r < 0, too, so that a{t) diverges in the "Big 
Rip" singularity at some finite moment of time in future. It is clear that, relaxing 
the requirement of scaling behaviour, it is possible to add some amount of dust-like 
matter to these solutions while still keeping the phantom behaviour of DE. However, 
the amount of 'fantomness' exhibited by them, i.e. the modulus of the minimal 
possible value of wde + 1 = 2/3g, is very small for small values of a (equivalently 
large ujbd), viz. 

a 

WDE + l>-- , (41) 

where the equality is achieved for n = 3. Thus, the conclusion is that polynomial 
potentials and scaling solutions in viable scalar-tensor DE models can lead to violation 
of the WEC, however, by the small amount ~ cu^^ only. 



5 Small z expansion and Solar System gravity tests 

Of course, scaling solutions considered in the previous section are very specific ones. 
Let us now study generic solutions describing DE in the scalar-tensor gravity. Since, 
as was explained in the Introduction, if DE crosses the phantom boundary at all, it 
has done it in a very recent epoch, at small z, it is natural to study the expansion 
of a generic solution in powers of redshift z. For each solution H{z), ^{z), the basic 
microscopic functions -F($) and U{^) can be expressed as functions of z and expanded 
into Taylor series in z: 



Fiz) 



1 + Fi Z + F2 z^ + ...>0 , (42) 



F[q 



Fo 

: ^luflU = ^u,o + Ui z + U2 z^ + ... . (43) 



Note that this expansion produces two parameters in each order of z which are inde- 
pendent of initial conditions and can be expressed through derivatives of -F($) and 
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[/($) with respect to $. The corresponding expansion for is: 



(Fo)-^$'(^) 



+ 2 + $'2 z' + ... 

A + ^ 6(Fi - F2 + %o - 1) - K^rafi + Ml) 



+ 




(4Fi - 2F2 + 6^]c7,o + 3r]„,o - 6) ■ z + . . . (44) 



with = 6 {VlDEfi — ^ufl — Fi). As we will see below, A^ > 0. In principle, the 
expansion (jUj) can be inverted to get 2($). From fl42|43p . all other expansions can 
be derived: 



which can be used in order to constrain parameters of our model. 

There are two types of observational constraints at small redshifts. The first of 
them includes those ones that follow from Solar System and other tests of possible 
deviations of the scalar-tensor gravity from GR at the present moment (z = 0); in 
particular, l^fo] < 10^'^ at the 2a confidence level from the last of Eqs. f|T3|) . Other 
constraints follow from cosmological tests and refer to the whole range of redshifts up 
to 2 ~ 1 and higher redshifts, depending on the nature of the test. In particular, if 
we assume that the present supernova data admit (or even favour) wde{z) < —1 for 
z < 0.3 as was argued in the Introduction, then ^ = ei + 2 €2 z + 3 63 z"^ + .... < 
for z < 0.3. In particular, we must have ei < that is equivalent to Eq. (I25l) taken 
at z = 0. 




1 + hi z + h2 z'^ + ... , 
1 + ei z + 62 z^ + ... > 

Wo + Wi Z + W2 z'^ + ... , 



go + gi z + 92 + 



(45) 
(46) 
(47) 



(48) 
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It follows from the substitution of these expansions into Eqs. f lTfS]) that 
1 



hi 
h2 



^{6- Sn^^o - Qnu,o - 4Fi + 2F2) 



3 

2^ 



+F| - 3F2^][/,o - ^i^2^^m,o + 4F2 + F3 - 5^]t/,o - Ml - 4fi„,o + 5 



1 + Wo 



Wi 



|Fi — 2F2 — 6{QdE,0 — ^Ufi) + f-^l ^DEfl 



3^DE,o{^ 



.(49) 
(50) 



3fi 



1 + 6 



Fi _ 1 
2 



4Fi - 2F2 + GQufl - 3^DE,o - 3 - 9 



mO 



n 



DE,0 



+ [If^ - 3F1F2 - ^ + 4Fl^]^,o + ^ + y Fi^^^o - 7Fi 



+ Fi - mnufi - -F2n^o + 4F2 + F3 - 5^] 



Ufl 



Ui + AVlDEfi + 1 



'iOS f^i _ n2 



4Fi - 2F2 + enu.o - 3f^z5i,.o - 3 



(51) 



The quantities Fi, F2, f2£i£;^o~^c/,o satisfy important constraints. For cu^^j^o; we have 
the expression 

_ 6{nDE,o - ^u,o - Fi) _ 

^BDfl - p2 - y^-^) 

which should be very large and positive, see (fT4|) . Therefore, we must have |Fi| -C 1 
and A2 ^ 6{Qde,o - ^Ufi) > 00 Moreover, since < GQdeo < 5 for positive U, 



\Fi\ < 



1/2 



< 10- 



(53) 



Thus, two cases are possible. In the first case, the further coefficients F2, F3, etc. 
in the expansion P2|) are all of the order of Fi, i.e. they satisfy the inequality fl53|) . 
too. In this case, the first derivative of F with respect to 2; or $ at the present moment 
is not atypical compared to other derivatives. Then, however, a possible amount of 
phantomness (the quantity defined in the end of the previous section) is also of the 
order of |Fi|, i.e. less than 1%. Such DE will be practically indistinguishable from a 
cosmological constant. 

Another possibility which admits 'significant phantomness', namely, m.m{l+ W£,e{z)) < 
—0.01, takes place if IF2I, IF3I, {^de.q ~ ^(/,o) a^^d so on are significantly larger than 



^Note that this condition is not satisfied in the recent paper ,67] that resuhs in unphysical nature 
of its best-fit solution for z < 0.2. 
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It is a matter of taste if one considers the second possibility as 'fine tuned' with 
respect to the first one; observations should finally tell us if significant phantomness 
does exist or not. In any case, it is clear that if we are interested in any prediction for 
a significant deviation from the cosmological constant variant of DE, we may neglect 
Fi as compared to F2 and other parameters (but not in those expressions where it 
enters multiplier). 

For <^ 1, all the expansions above simplify significantly, and we have, in 
particular. 



1+Wo 



2F2 + 6{nDEfi - nu,o) 



From fl54|) . the necessary condition to have phantom DE today reads 

'd^F\ Fo 



(54) 



d$2 



3 {n 



DEfl 



-n 



< -1 



u,o) 



(55) 



In particular, F2 < 0, because ^de,o — ^u,o > as discussed abovej^ In the same 
limit, hi < 0, i.e. the WEC is violated for the total effective energy-momentum tensor 
for matter + DE (for the whole right-hand side of Eq. (fTS!) ). if the stronger inequality 
is fulfilled: 

F2<~il + nnEfl-2nu,o) . (56) 
The expression for wi becomes in the limit ^ 1 
1 



Wi 



+ 



'^^DE,0 
2 

^DE,0 
1 



1 + 6-^ ) {2F2 + 3 + 3nnE,o - 6nu,o) - 9-^ 



Ufl 



-^m,0 + 4 



5Qu,0 + ^^DEfi + 1 + 



Ml 



{2F2 + 3nDE,0 + 3 - Qnufl) 
'^^^DEfi 

9{wo - ifriDEfl -8-6w'^ + 9Q 



(57) 



m,0 



ISQufi + Wo (26 - 9n„,o + I8n 



U,0 



+ 



[-Qwo + lAQufi + 31)^,0 (-3^^(7,0 + 2u7o)) 



n 



DEfl 



n 



DEfl 



(F3 - Ui 



{51 



where the last equality (l58l) is obtained using (15^ . Since F3 and ui are free parameters 
determined by a concrete choice of F{^) and U{^), it is well possible to have wi > 
for Wo < —1 in order to realize a smooth phantom divide crossing at small z for DE 
in the scalar-tensor gravity. 



^In this place our considerations intersect with those of the recent paper [68] . The condition ((55)) 
corresponds to /3o > 1 in the notation of that paper. The other case mentioned there, f3o ^ —1, does 
not lead to significant fantomncss of DE with our definition of the DE energy- momentum tensor. 
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9o 



We obtain further 
= FAl- 



1 



X 



- 4F2)(y - 1) - (5Fi - 2F2 + 6fi^,o - 3fi^ij,o - 5)Fi 

2 



■ (59) 



Expression (l59l) simplifies considerably for ^ 1: 

F2 



90^ Fill 



Fill- 



(60) 



Note that go and Fi have the same sign for the case of (significantly) fantom DE. 
Finally, for the post-Newtonian parameter /3p7v and 7PAr we have 



iPN 



1 + 
1 



BD 



4(2cubz5 + 3)(cubb + 2)2F' 
1 



ubd + 2 

Taking the zeroth order in z, we obtain : 



Ppn,o 



1 + 
1 



^BD,1 



1 



4{2uJBD,0 + 3){uJBD,0 + 2yFi 

1 



For ^ 1 and \FiUi/A'^\ ^ 1, we have: 



(61) 
(62) 

(63) 

(64) 



UJpDfi — Q{^DE,0 - ^U,o) ^2 



1 

F? 



J- 2 



(65) 
(66) 



which finally yields 

lPN,0 = 1 



F? 



6{n 



DE,0 



-n 



/3. 



PN,0 



F?F, 



Ufl) 



72 (fi 



DE,0 



(67) 



Now we can use f l50|52ll60|l67|) to extract information from the Solar System con- 
straints. First we note that all the Solar System constraints f lT^ . flT^ are satisfied 
for a sufficiently small (less than the upper bound flS5]) ) independently of any 
requirement concerning the present DE equation of state, in particular, whether one 



16 



has today DE of the phantom type or not. From (1521) . it is seen that 7pAr,o ~ 1 and 
7PAr < 1 for any viable scalar-tensor model of DE having ujbd,o positive and large. 
For significantly phantom DE, other constraints follow from ([6] 



Ppn,o > 1 , (68) 
- = 6 ilDEfl — 4. (69) 



PpN,0 — 1 ^2 

Therefore, for significantly phantom DE 



-4<|^^^<0. (70) 



PN,0 



It is possible to invert formulas ( l60ll671) and express Fi, F2 and flDE,o — ^u,o 
through the Solar System observables 7piv,o, Ppn,o and go: 



Fi = 


90 

7 - 


F2 = 


-2 g'o ■ 




1 2 


l + Wo = 


1 2 
-3^0 



7-1 



l-4(/3-l) 
P-1 



[7-l-4(/3-l)]2 

7-1 
[7-l-4(/3-l)]2 

4(/?-l) + 7-l 
^DEfl [7 - 1 - 4(/5 - 1)]^ 



(71) 
(72) 
(73) 
(74) 



Thus, in principle, it is possible to test wq < —I in the Solar System as was recently 
discussed in [68] • However, this may be very difficult to perform in practice since 
the small parameter Fi enters quadratically into 7pAr,o ~ 1 and f3pN,o — 1 while not 
appearing (in some limit) in 1 + wq. So, a rather significant fantomness of DE is 
typically accompanied by very small deviations of the post-Newtonian parameters 
from their GR values. E.g., let us take floEfl = 0.7, flufl = 0.6, 7p7v,o"~l = — 2-10~^ 
(the latter being marginally possible at the 2a level) and Wq = —1.2. Then, from Eqs. 
( 16711541) . we get |Fi| = 3.5■10-^ F2 = -0.51 and /5pjv,o - 1 = 0.85 ■ 10"^ - an order of 
magnitude below the present upper limitJzl Further, wi = 2.9 (F3-M1) -0.70. Finally, 
in this case \go\ = 0.95 • 10~^ that is an order of magnitude larger than the 2a upper 
limit following from the last of Eqs. (fT3|) ! So, if this upper bound will be confirmed, 

has to be decreased by an order of magnitude which results in Ppn,o — 1 and 
7pa',o ~ 1 being on the level of 10^''. For comparison, in the extreme opposite case 

= with the same values of floEfl, 1pn,o and Wq, we get = 9.2 • 10"'^, F2 = 
-2.3, /3piv,o-l = 0.55•10-^ Wi = 2.9 (F3 -mi) + 1.5 and even larger |^o| = 1.9-10-2. 

This shows also that the measurement of Gefr,o/G^efr,o is the most critical among 
Solar System tests of scalar-tensor DE since this quantity is proportional to the first 
power of the small parameter (apart from the exceptional case {(PF / d^'^)^ = 1 



^Let us emphasize once more that the formulas (|60l67p were obtained under the assumption 
1^2! ^ l^il- For this reason, they are not appUcable, e.g., to the scahng solution (|40l) for which 
F2 = —Fi/2 > 0, 7PAr — 1 — ia in the hmit a <^ 1, and f3pN = 1- 
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which does not lead to the present phantom behaviour of DE). Also, to determine Wi 
which is necessary in order to consider the possibility of phantom boundary crossing, 
the determination of {(f\nGes/dt'^)Q is required, something that is hardly possible. 
Thus, testing the phantom behaviour of scalar-tensor DE in the Solar System may 
be much more difficult than in cosmology. 

6 General reconstruction of F{z) 

We consider now the reconstruction of F{z) for given Hubble parameter H{z) and 
potential U. As can be seen from fl24f2ip . H{z) is a functional of w{z). We will 
consider several cases corresponding to phantom divide crossing at very small redshifts 
as favoured by the latest observations. It is desirable to derive some general properties 
of the behaviour of our system before embarking on the study of specific models. 

It is easy to derive the following general property shared by a large class of models: 
if we require that 0'^ > and also U > 0, then the following inequality must be 
satisfied 

F' /(I -\- z)^ F' \ F 

[^^y -^^ + '^)+Y^>^^ + ■ (75) 

Note that the r.h.s. of (1751) tends to one in any model for which wde < ^ and it 
will tend to this asymptotic value quickly for wde < —0.5. We can consider several 
particular cases: 

• It is seen from f l75|) that models which at some redshift ^ 1 satisfy F'{zm) = 

and < F{zm) < 1 sufficiently small, must also have (p''^{zm) < 0, such an 
example is actually illustrated in Figure 1. 

• When F —>■ for z z^ while at the same time F'{zm) 0, it follows from 
the expression ffTOj) for (f)'^ and using ([8]) that 0'^ ^ oo. 

• Finally, an interesting case is provided when F{zm) = F'{zm) = for some z^- 
We see first from ([8]) that — > $'^(2^), where ^'"^{zm) is a (small) negative 
number. Let us assume that ^ is bounded when z Zm- In that case the 
inequality ( 1751) cannot be satisfied, hence ^'^(^m) < (allways assuming f/ > 0) 
and actually we have 0'^ — > —00 for be checked directly with 
the definition of 0'^. As a consequence, if 0'^ > as 2; — Zm then we must have 

1 -^r I — >■ 00 and if in addition is bounded, then 0'^ — ^ 00. 

Some of these properties are illustrated with Figures 1,2,3. 

Let us consider now the space-time background evolution encoded in the quantity 
h{z). From a theoretical point of view, as is seen from (12^ . a given function ^(2;) 
implies a corresponding functional form h{z), and conversely using (1221) . Future data 
are expected to mesure Dl{z) = Hq^cIl^z), and therefore h{z) = Hq^ H{z), with 
high precision. In the meantime, we can try some particular expressions w{z), or h{z), 
parametrized with the help of a limited number of free parameters, an attitude which 
turns out to be very fruitful. We will consider both constant and variable equation 
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of state parameter w. Observations suggest that w can be varying and we will model 
this variation using the following two-dimensional parametrization of the equation of 
state parameter suggested in [17], [T8] ) 

w{z) = {-1 + a) + p (l- x) =wo + wi . (76) 

X ~|~ z 

where x = We will sometimes compare our results with the parametrization 
suggested in [19j 

e{z) = Ao + A^{l + z) + A2 (1 + zf . (77) 
By definition, e(0) = 1, hence Aq + Ai + A2 = 1. Note that we have for (j76ll 

e{z) = (1 + 2)3("+/3)e-=^^TTi . (78) 

Actually, rather than delimiting some restricted bounded domain in the parameter 
space, some observations on small redshifts z < 0.35 single out a prefered direction. 
Variations along this direction are essentially unconstrained while variations normal 
to it are most efficiently constrained. This is the case for baryon oscillations data 
which constrain h{z) and therefore w{z) as follows 



1/3 




dz ' 




.^1 


/o Kz)\ 



2/3 



< 0.469 ±0.017 , (79) 



with zi = 0.35, flrn,o = 0.3. 

When we use the parametrization (1761) . eq. (!79|) translates into the constraint at 
the l-cr level in the parameter plane a, (3 (equivalently wq, wi), 

a + 0.112 (3 = 1 + WDEfl + 0.112 wi = 0.23 ± 0.20 . (80) 

As said in the Introduction, we see that these constraints allow for a Wq close to (and 
possibly slightly lower than) —1 if we take a constant equation of state, while a Wq 
significantly lower than —1 is allowed but it requires a steeply increasing w near z = 
in agreement with recent analysis of the data. It is natural to use the constraint ( |80l) 
when we consider the small z expansion of our scalar-tensor model quantities and of 
our fit (EE]). 



7 Reconstruction for a constant potential U 

Setting U = 0, we consider either wde < constant, or else wde{z) of the form 
( I76ll or (1771) . whose variation is parameterized by two parameters. This extends the 
investigations in |49j, where only the case wde = —1 was considered. The interesting 
result obtained there was that F vanishes, and hence the theory is singular, at very 
low redshifts z ~ 0.66 for -F'(O) = and at only slightly higher redshift z f» 0.68 if one 
fully exploits the possible initial conditions allowed by the solar system constraints. It 
is therefore interesting to investigate how this result is affected when we take different 
equations of state. In particular we will also consider here equations of state with 
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Wde < — 1- As SNIa observations point to a varying equation of state wliicli is of 
tfie pliantom type on very low redsliifts z < 0.3, this case is considered too. For 
the polynomial expression (1771) . for definiteness we will use the best fit to the "Gold 
data set" consisting of 156 supernova which was proposed in [2^ (with the priors 
^m,o = 0.3 and VL^ q = 0): 

= -5.94 ±3.61 ^2 = 2.39 ± 1.47 . (81) 

These values are similar to those earlier obtained in [20] for the same dataset. As we 
can see, (|5T]) has large uncertainties. The corresponding equation of state is today 
of the phantom type but rapidly increases to a positive value wde — 0.1 at z ~ 0.8 
and then decreases to its asymptotic value wde = — | but it is still slightly positive 
at z = 2. We first note from (|5^ that F2 satisfies 

F2 = ^ {wDEfl - 1) ^DE,0 + 3 flufl . (82) 

Therefore F2 < for ?7 = as long as wdea < 1 while F2 decreases with decreasing 
wde,o- The results obtained numerically are shown in Figure 4. It is seen that we 
keep essentially the same picture, in particular F vanishes more rapidly when W£,e 
decreases, hence the problem is even more severe for phantom DE with constant wde- 
Even the fit (IHTj) where wde varies substantially at low redshifts yields essentially the 
same behaviour. The same result as for flHTj) is obtained using the parametrization 
fl76|) with a = —0.377 and j3 = 2. This is similar to earlier results showing that a large 
variation of the equation of state around wde = —1 on low redshifts can result in 
essentially the same dL{z) pn[33j for z ^1 and one can understand from the general 
expression eq. (!3T|) why all cases displayed in Figure 1. will have basically the same 
behaviour regarding the evolution of F{z). The initial condition Fq = Fi = and 
$'(0) 7^ means physically that the theory corresponds strictly to General Relativity 
today, i^^Do ~ 0. The solar system constraints allow for a very small nonvanishing 
Fi, see eq. fll5p . and the corresponding change in Zmax is very marginal as can be seen 
on Figure 4. To summarize, looking only at z = one could think that models with 
u!de,o = —1 and U = are allowed, however we see that the cosmological evolution 
of such models leads to a singularity at very low redshifts generalizing the results 
obtained earlier [19]. It is clear from these results that a scalar-tensor theory of 
gravity with vanishing potential U is definitely excluded by the Supernovae data. 

The next natural step is to consider constant (nonvanishing) potentials U. Of 
course, one does not expect such a behaviour to be relevant up to very large redshifts, 
but it is certainly a sensible approximation to start with on small redshifts z < 1.5. 
We note from (!82|) that F2 < for Wde,o < By taking Quq ^ ^de,o one has 
the smallest possible F2 for given weie,o and these are the cases which are found to 
have the largest range of validity. We find that many models are allowed which are 
perfectly viable on small redshifts z < 1.5. Typically, these models become singular 
at some higher redshift well beyond z = 2. There are several possibilities: either it is 
the quantity F that is vanishing first or it is the quantity 0'^ that vanishes first. As 
discussed in Section 6, when F = and F' 7^ 0, the quantity 0'^ diverges, 0'^ 00. 
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Models for which both F and F' vanish together can be considered as a hmiting case 
which gives the largest possible redshift. If we change very slightly the potential U 
so that when F' = 0, F > 0, one gets 0'^ = earlier. These different possibilities are 
displayed in Figures 2,3. We have checked the behaviour of such models with wde 
starting below —1 as favoured by SNIa data and with wde quickly becoming larger 
than -1. To summarize, we find that models with f2[/o — ^ 0.7 and an equation of 
state with W£,e —1 will become singular at arbitrarily high redshifts. We should 
remember that Qu^ = 0.7 together with wde = —1 gives back General Relativity 
[F = 1). In all other cases some maximal redshift is found where the model becomes 
singular. 

In all our numerical calculations, we neglected radiation since its energy density 
is very small at redshifts of interest. However, even in principle its presence cannot 
prevent the occurence of the singularity at the moment when Ge// changes sign (0'^ = 
in our case) whose generic (anisotropic) structure is independent of the matter 
equation of state (see [55] in this respect). 

8 Asymptotic stability 

The next step is to consider nonconstant potentials. As we have shown very generally 
in Section 5, such models are consistent with DE of the phantom type today and 
phantom divide crossing at small redshifts. An example of the reconstruction of a 
model with phantom DE today and phantom divide crossing at 2; ~ 0.3 was presented 
in [50j. However, in this paper the reconstruction was implemented only for small 
redshifts z <2. As was shown in the previous section, one has to consider the large- 2; 
behaviour of a model, too, to prove its viability. 

In contrast to the case with constant or vanishing potential f/($), a growing 
potential U{^) allows for the construction of scalar-tensor DE models which are 
viable for all redshifts and evolves according to the fit (1761) 

w{z) = const = wo + wi = -l + {a + f3) > -1 z > 1 . (83) 

As a particular example of such a model, let us assume that the scalar-tensor gravity 
approaches GR sufficiently fast at the matter-dominated stage and that DE 'tracks' 
matter: 

F ^ F^ = const, |F| < HF^, \F\ < H^F^, (84) 
(X {1 + zf K z < Zeg , (85) 

where Zeq is the redshift at the matter-radiation equality. In this way we recover in 

2 

the past the usual behaviour a ccts. However, the constant value F^o < Fq may not 
be too small compared to Fq. In order to satisfy the BBN constraints [691 EO], the 
following inequality is required 

(Fo - F^)/Fo < 0.1 , (86) 
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which can be easily satisfied. Indeed, as can be seen from (1751) . we have for constant 
F and [/ > 0, 0'^ > 

Fo-F^< — ^^Fo . (87) 

Pm + PDE 

On the other hand, the assumption Fq > F^o provides a good matching to the small-z 
expansion derived in Sec. 5 with F2 < and a very small Fi. In view of (1861) . this 
matching should occur at a sufficiently small z, too. 

This scaling behaviour of DE corresponds to the asymptotic solution (l85l) with 
wq = —wi, equivalently a + j3 = 1. As is well known, it can be obtained by taking 
an exponential potential f/($) at $ — ^ 00, i.e. at large redshifts: 

U oc eVwl^* ^ (88) 

where flu,oo is the constant asymtotic value of flu at z ^ 1 which is a free parameter 
formally. However, actually it should be small (less than a few percents) to obtain the 
correct value of the growth factor of density perturbations during the total matter- 
dominated stage. The total DE energy density in terms of the critical one 3FqH^ 
is 

nDE,oo = 21^C/,oo + < 1 . (89) 

-TO 

Note that the term in (15^ always has the same equation of state as the main 

background matter. Thus, to obtain behaviour different from wq = —Wi at large 
z during the matter dominated stage is possible if Fq = F^o only which requires 
additional fine tuning and is not natural. 



9 Conclusions 

In this work we have considered the viability of scalar-tensor models of Dark Energy. 
We have used different types of observations: Solar System constraints which con- 
strain the model today, and other data like the supernova data which constrain its 
cosmological evolution, in particular the time evolution of the DE equation of state 
parameter wde- We were interested specifically in models which violate the weak 
energy condition on small redshifts z < 0.3 and in any case today. We have found the 
formal general integral solution for F{z) when we reconstruct it for given H{z), which 
can be obtained from the d^lz) data, and for given U{z). This general solution allows 
immediately for an integral representation of scaling solutions. We have constructed 
scaling solutions and shown that they exist in models with F = with a =constant, 
in these models the Brans-Dicke parameter ljbd is constant. Only for nonzero po- 
tentials can these models have scaling solutions with constant wde < —1- However, 
it is shown that for these models jw^i^; + 1| is very small with |1 + ~ ^^bd- 

We have further performed systematically the small z expansion of the theory and 
used it to extract various observational constraints. We recover that a large positive 
^BD,o requires |Fi| <^ 1 where Fi is the first derivative today of We find that a 
significantly phantom DE today {wq < —1.01) implies that F2, the second derivative 
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today of must be negative and not small, thus, significantly larger than Fi by 
modulus, fiowever, while necessary this is not a sufficient condition for phantom DE 
today, for example for vanishing potential F2 < whenever wdea < 1- 

The condition IF2I ~ 1 while the Solar System data require < 10~^ (i.e. 
anomalously small) is the only 'fine tuning' required to get a significantly phantom 
DE at the present time in scalar-tensor gravity. Note that, since the derivatives Fi{z) 
are not parameters of the effective microscopic Lagrangian (1) but depend on initial 
conditions in the early Universe too, it could be even better to call this a "cosmic 
coincidence". Our point of view is that the condition IF2I ~ 1 cannot be excluded by 
pure thought, only observations will possibly do it: in the absence of this condition 
when all Fi are of the same order as Fi, the general prediction is that the amount of 
possible phantomness in scalar-tensor DE models is very small, less than 1%. 

As for the solar system constraints, we have shown that the Post-Newtonian pa- 
rameter 7pAr,o must satisfy 7pAr,o — 1 — —'-^bda < 0' ^^^^ ^ general requirement 
for the viability of our model. On the other hand, a significantly phantom DE to- 
day implies PpNfi > 1- Combining those results, we find that the negative quantity 
does not depend on Fi and can be expressed in function of F2 and 1 + WDEfl- 
This would enable us to find F2 from the Solar System constraints provided we know 
1 + wde,o from other cosmo logical data. Still for phantom DE today, we find that 
Geff,o has the same sign as Fi. So, a measurement of the sign of Gefr,o would give us 
the sign of Fi if we know we have 1 +WDEfl < 0. On the other hand, a measurement 
of both quantities with opposite signs would rule out phantom DE today though 
measuring the sign of Fi can be hard to determine observationally. However, due 
to the already confirmed smallness of Fi, connection between cosmological and Solar 
System tests of dark energy is rather one-way since it requires much greater accuracy 
from the latter ones for the determination whether DE is phantom at present or not. 
While a positive detection of /^pat^o > 1 is a strong argument for the phantom DE 
at present, a negative result (no measurable deviation of Ppn,o from unity) tells us 
nothing regarding DE properties. Also, the Solar System tests are clearly unable to 
provide Wi in any reasonable future, so no information about the possibility of the 
phantom divide crossing may be expected from them. 

We have also considered numerically the reconstruction of F for various wde, 
including constant as well as varying equations of state of the phantom type. Gener- 
alizing results obtained in [19] for a pure cosmological constant, we find that models 
with a vanishing potential (see Figure 3) are ruled out and lead to a singular be- 
haviour for z < 0.66. While it is clearly possible to have phantom DE today with 
U = without conflicting with the data, the cosmological evolution of these mod- 
els rules them out. For models with constant nonvanishing potentials, which can be 
considered as a good approximation on small redshifts for more general models with 
varying potentials, we find that it is easy to have models in agreement with observa- 
tions on small redshifts z <2. However, it is interesting that these models generically 
have a maximal redshift where they become singular. So, to construct a scalar-tensor 
DE model having a sufficiently long matter-dominated stage, a non-constant poten- 
tial U{^) is required, and we have presented an example of such a model where DE 
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tracks matter at large redshifts. 

Therefore, the final conclusion is that the generic scalar-tensor gravity ([1]) with 
the two functions F{^) and derived from some underlying theory (e.g., from 

brane cosmology) or taken from observational data, has enough power to provide 
internally consistent cosmo logical models with a temporarily phantom DE (at present, 
in particular) and with a regular phantom divide crossing in the course of the evolution 
of the Universe. 
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5 10 15 20 25 30 

z 

Figure 1: Several quantities are displayed for the model with parametrization (176!) 
with a = 1 + Wo = -0.2 and /3 = Wi = 0.4, while Fi = and = Tt^ = 0.97585. 
The curves shown represent the following (rescaled) quantities in function of redshift 
z: (sohd), 10 x (dotted), 10 x 0'^ (dashed) and 0.1 x uj^d (dot-dashed). It is 
seen that lobd and become negative at z ~ 10. The model remains valid beyond 
z ~ 10 untill z ~ 18 as long as ujbd > — f or equivalently 0'^ > 0. So, we see that the 
model remains valid for a large interval where < 0. Of course, it is impossible to 
reconstruct $ in this interval using the Z = 1 parametrization. Note that 0'^ becomes 
negative before Zm ~ 22 where F'{zm) = 0, in accordance with (1751) . 
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Figure 2: The quantity is displayed for the parametrization (175]1 with a = 
1 + Wq = —0.2 and P = wi = 0.4 and Fi = 0. We have the following values for 
= ^ from bottom to top: 0.9758, 0.975824492, 0.975826, 0.97587. The 
second curve has its minimum at F = and is superimposed on the third curve 
which has its minimum at F = 2.4 x 10^'^. 
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Figure 3: The quantity is shown for the same models as Figure 2. The short, resp. 
long, dashed curve corresponds to = w = 0.975824492, resp 0.975826. 
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Figure 4: The function is shown for several models with vanishing potential f/ = 0. 
The solid lines correspond to the initial conditions Fq = while the dashed lines 
correspond to the maximal values allowed by the solar system constraint ujbd,o > 4 x 
10^. We have from left to right the following equation of state parameter w: —2, —1.5, 
polynomial expression flHTl) . —1, —0.5. It is seen that the limit of regularity of these 
models corresponds to very low redshifts, 0.566 < z < 0.663 for —2 < w < —1. Note 
that the polynomial expression flHTj) represents crossing of the phantom divide. 
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APPENDIX 



We summarize in the table below, the integrability of scaling solutions using (|3T|) for several exponents n = 87. We use the 

following abbreviations: 

a.n.i.: analytically non integrable 

Elliptic, log, Hypergeo : solutions expressed in terms of resp., elliptic, log (or argtanh) and hypergeometric functions 
P^"'\x): polynomial in x of degree n 



n 
w 


< 

< -1 




-1 


1 

-2/3 


2 

-1/3 


3 



4 
1/3 


5 

2/3 


6 
1 


7, 8, 9 
2/3,5/3,2 


Ao 

Ai 
A, 


a.n.i. 
a.n.i. 
a.n.i. 


Elliptic 
log 
log 


Elliptic 
Elliptic 
Elliptic 


log 
log 
log 


2 

y^(l+k)x 
2 


k+l/x 


Elliptic 
Elliptic 
Elliptic 


Elliptic 

2 l+kx''^ 


Hypergeo 
Hypergeo 
Hypergeo 


3A/(l+fc)x3 
2 


x^ k+l/x 
P(5){a;) 


3 Va;=^+fca;6 
2 il+kx^)(-l+2kx^ 




x''y/ k+l/x 


9 x'^^/x'^+kx*' 


A 


a.n.i. 
a.n.i. 


a.n.i. 
a.n.i. 


a.n.i. 
a.n.i. 


a.n.i. 
a.n.i. 


1 

3(l+fc)x2 
2 

45(1+A;)a;5 


1 _ Ak 
3x'^ 3x 
2 1 256fc** 
45x5 3142: 


a.n.i. 
a.n.i. 


a.n.i. 
a.n.i. 


a.n.i. 
a.n.i. 
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